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We reformulate the Tomonaga-Luttinger liquid theory for quasi-one-dimensional Fermion systems with many 
subbands across the Fermi energy. Our theory enables us to obtain a rigorous expression of the local density 
of states (LDOS) for general multichannel quantum wires, describing how the power-law anomalies of LDOS 
depend on inter- and intra-subbands couplings as well as the Fermi velocity of each band. The resulting formula 
for the exponents is valid in the case of both bulk contact and edge contact, and thus plays a fundamental role 
in the physical properties of multicomponent Tomonaga-Luttinger liquid systems. 
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I. INTRODUCTION 

Bosonization is one of the most powerful techniques for 
describing the properties of one-dimensional (ID) interacting 
electron systems. In ID systems, even a slight interaction be- 
tween electrons strongly affects the quantum nature, result- 
ing in the occurrence of Tomonaga-Luttinger liquid (TLL) 
statesji^ TLL states exhibit power-law anomalies in phys- 
ical quantities, as predicted by the bosonization theory.^ A 
prominent example is the power-law singularity of the single- 
particle density of states D{E, T) near the Fermi energy, E-p, 
represented by D{E,0) oc \E - Ep]^ and D{Ef,T) oc T^, 
with E and T being the energy and the temperature, respec- 
tively. The value of A, called the TLL exponent, is dependent 
on the interaction strength^ and other parameters character- 
izing the ID system.i§^-ii Recently, it has been suggested that 
a continuous variation in A can be produced by an external 
fieldr^— this implies artificial control of the transport prop- 
erties of quasi- ID conductors, since A governs the power-law 
behaviors of the differential tunneling conductance^ dl/dV cx 
at high bias voltages {eV ^ ksT) and the temperature- 
dependent conductance G{T) cx at low voltages [eV <C 
kBT). 

Experimental realizations of TLL states encompass 
various systems showing highly anisotropic conduc- 
tivity: metallicr^^ semiconductingr^^^^ and organic 
nanowires,-^22 and carbon nanotubes^iSill^ are a few 
examples. These actual quasi-lD conductors possess a finite 
cross-section, thus exhibiting a finite number of transmission 
channels in the transverse direction (except for a limited case 
in which E-p is small enough for only the lowest subband to 
be involved). The presence of multiple channels at Ep causes 
inter-subband scatterings. Furthermore, different channels 
can have different Fermi velocities, i.e., the slope of the 
dispersion curve at E-p (see Fig. [U, and thus contributions 
from each channel to the TLL exponent differ from each 
other Theories of multichannel TLL have been developed for 
the Hubbard model in the presence of an external magnetic 
field,——, where the discrepancies in the Fermi velocity 
between up- and down-spins are taken into account. A 
similar issue was also discussed in the study of quasi- ID 



Bose gasesi^ The effect of inter-subband scattering on the 
TLL exponent has been investigated in connection with the 
TLL behavior of multiwall carbon nanotubes.^- However, to 
the best of our knowledge, singular behavior in the density 
of states D{E, T) remains unresolved for multichannel TLL 
systems with the coexistence of inter-subband scatterings and 
Fermi-velocity variations. Hence, the rigorous expression of 
D{E,T) in multichannel TLL systems is desirable for de- 
scribing the transport properties and photo-emission spectra 
those will be experimentally observed in actual quasi- ID 
conductors. 



In this paper, we reformulate the multichannel TLL theory 
in order to derive the anomalous energy- and temperature- 
dependences of the local density of states of quasi- ID 
Fermion systems. Cases of locations both far from the bound- 
ary and close to it, which correspond to bulk contact and end 
contact of the transport properties, respectively, are discussed. 
We demonstrate clearly how the TLL exponents of multichan- 
nel systems depend on mutual interaction and Fermi veloci- 
ties. The resulting formula for the exponents, as well as the 
theoretical framework we have established, will provide clues 
to exploiting the effects of subband couplings and Fermi ve- 
locity variations on the nature of TLLs in real ID systems. 



The paper is organized as follows. In Section II, the mul- 
tichannel TLL theory is developed for A^-channel quasi- ID 
Fermionic systems with different Fermi velocities. The local 
densities of states far from and close to the boundary are cal- 
culated in Section III. As a simple example, in Section IV, the 
case of spinless Fermions is discussed and the A^-dependences 
of the exponents are obtained for long-range interaction lim- 
its. The paper closes with a summary in Section V. In the fol- 
lowing, the unit ft = fcs = 1 is used, unless explicitly stated 
otherwise. 
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FIG. 1: Sketch of the energy dispersion of the present system, with 
energy bands cross the Fermi energy, Ep. The Fermi velocity and 
the Fermi wavenumber of the i/-lh band (v = 1, - ■ ■ , A'^) are denoted 
by vpi, and fcp,y, respectively. The symbol p = + (— ) indicates a 
one-particle state moving toward the right (left). 



II. MULTICOMPONENT TOMONAGA-LUTTINGER 
LIQUIDS 

A. Bosonization 

We consider a quasi- 1 D Fertnion system with iV ID energy 
bands cross the Fermi energy, Ep. The band structure close to 
E-p is schematically shown in Fig. [T] Here, the Fermi velocity 
and the Fermi wavenumber of the i/-th band {ly — 1, • • • , N) 
are denoted by vpi, and kp^, respectively, and the one-particle 
state moving to the right (left) are indicated hy p = + (— ). 
The kinetic energy of the Hamiltonian, "Hk, is expressed by 



N 



,y=lp=± k 



(1) 



where the one-particle energy and the wavenumber k are mea- 
sured from Ep and pfepi^, respectively. In Eq. cl ^ ^ de- 
notes the creation operator of the Fermion with wavenumber 
k, branch p, and band index i/. 

Let us introduce the density operator of the p-branch of the 
ly-th band, defined as 



= T,k 4,p,^Ck,p,t. : • • • q = 



, (2) 



which satisfies the commutation relation 

[Pp,i^{-^q),Pp'y{q')] = 5pp>5yi,>5qq>pqL /{2t:). In terms 
of 

Pp,i^iq)' is expressed by^ 



N 



^k = X X PpAq)ppA-q)^ 



(3) 



where L is the length of the system. The most general form 
of the mutual interaction between Fermions leading to the iV- 
component TLL is written as^ 



jy,l/'=l p.q 

+ 9i{v-, iy')pp,uiq)Pp,u'i-q)y 



(4) 



The matrix elements 52 (i^, J^') and 54 (z^, ly') depend on the de- 
tails of the model we consider. Specifically, the case with 
92 = 94 corresponds to the model for multiwall carbon nan- 
otubes considered in Ref. [s^ As an example, we will discuss 
the case of the spinless Fermion in SectionlTVl 

We introduce the phase variables 9^{x) and (j)^{x) {v = 
1, • • • , N), defined as 



2Trpx 
L 



Np,u 



— z_.p—^PpAq)~P2^-P''' 

2TTpx 



(t)u{x) 



(5) 



g#0 



-Np,y 



2m e "^-^ , , TT I 
— 2^P—;;-PpAq)^P-;^^-P,y (■• (6) 

' 97^0 " 



where [Qp.,j,Npi y] = \5ppi5i,^i. In the summation in terms 
of q, the ultraviolet cutoff exp(— Q!|(7|/2) is implicitly in- 
cluded. The phase variables satisfy the commutation relation, 
[6^{x),(t)^,{x')\ = i2'K5^^,e{x - a;') for L ^ 00, with 6{x) 
being the conventional step function. In terms of the above 
phase variables, the Hamiltonian is written as 

1 " r 

(7) 

where = —dx<pi^/{2n). This is the general form of the 
phase Hamiltonian expressing the A^-component TLL. The 
symmetric matrices K and V are defined as follows: 

{K~^)^,y = 2TrvFuS„„' + gi{v, v') ~ 92(^1 ^'), (8) 

V„u' = -^o^u' H . (9) 

2tt 47r^ 

The Fermion operator, defined by ipp.^{x) = 
VL^i^e^'^^Ck^p^^, is related to the phase variables 

as 



'^P,>^i^) = /T^ cxp <^ i— = [et,{x) +p(l)^{x)] > , (10) 
y/ZTTa I v2 ) 

where rj^ expresses the Majorana Fermion operator satisfying 

T]„ = nil and {-qu.-nv'} = 25^,,,. 



B. Diagonalization 

The Hamiltonian given by Eq. Q has a bilinear form with 
respect to dx(t>i> and dxOv, and thus can be diagonalized by the 
standard unitary transformation, as shown below. 
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The equations of motion of the phase variables derived from 
Eq. Q read as 



^^j^ T, x). The contribution from the i^-th band, 

Di,{uj, T, x), is given by 



9 jj_y Q 

dt dx^ ' 



d_ 

dt 



(11) 

(12) 



with6» = {91,62,- ■■ ,0jv)^ andn = (ni,n2,-- - ,UnV- 
Here, the energy eigenvalue, w = v\k\, and the eigenvector X 
corresponding to it are determined by 



{v^K - V)X = 0, 



(13) 



whose solutions are denoted by vj and Xj (j = 1, 2, • • • , N). 
The eigenvectors are normalized as {Xi, KXj) ~ Sij. 

To obtain a concise representation of H, we define the uni- 
tary transformation as 



e = x&, (14) 

n = KXa, (15) 

(ei,e2,--- ,eAr)T and h = 



where = 
(Si,S2,-- - The N X N matrix X consists of 

the set of eigenvectors Xj as X = {Xi,X2,--- ,Xn), 
and satisfies X"^KX = 1. Under the transformation, 
['dj{x),'E.jr{x')] = i5jjr6{x — x'). By using the new 
variables, we obtain an alternative form of H, given by 



N 

and that of the field operator defined in Eq. (fTUt . 



(16) 



V2t 



»} 



(17) 



where = ^f^i(i^X),,$j 



III. DENSITY OF STATES 

In this section, we discuss the local density of states 
D{uj,T, x) with u! = E — Ep for lu ^ Ep, where x denotes 
the position along the ID direction. As noted earlier, the TLL 
exponent that characterizes the singularity of D{uj, T, x) near 
Ep is X dependent. From a practical view, it is specifically 
interesting to study the semi-infinite system with its end at 
the origin^** and discuss the cases with x — > and x oo, 
which correspond to the end contact and the bulk contact, re- 
spectively. In the following argument, we therefore derive the 
TLL exponent for both cases, as well as the explicit forms of 
D{uj, T, x) as functions of u and T. 

The local density of states is given by the summa- 
tion of the contribution from each band: D{uj,T,x) — 



1 1"°° 

D,{uj,T,x) = — tc'"*({Vi(x,0),V.(a:,i)}), (18) 



where tjjuix.t) 



Here, {A, B} = AB + BA and (• • • ) means the thermal av- 
erage. The quantity in the integrand in Eq. (fTsT l for x ^ cx) is 
proven to be 



{{iplix oo,0),i:^{x -> oo,t)}) 
Tra L J 



(19) 
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N 

n- 

j=i (1 - ivjt/a) "-J 

(20) 



(21) 



where the superscript (6) means the case of a "bulk" contact. 
Similarly, the counterpart for the "edge" contact case a; 0, 



labeled by (e), obeys Eqs. (fT9]l-(l2ni with Y^^j replaced by 



=2nv, [{KX),J 



(22) 



The derivations of Eqs. (fT9]l-(l22]i are shown in Appendix A. 
We are ready to obtain the rigorous expression for 

D^^''\i^,T) = D,{ij,T,x ^ oo/0)4 Since ^''^"^(0,0) 

vanishes as long as J2f=i ^u^j"^^ > 1' we subtract it from the 
result and take the limit a — > 0. Eventually, we attain the 
desired formulas: 



1 



t< e' 



iujt 



f TTTt \ ^^ = ' -'^ 



27r^a J-ao I \sinh7rTi^ 

N , N 



X 



n 



1 



n 



- 1 



1 




Equation ( |23] | implies that the uj- and T-dependences of 
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D^y'^'^^ (u, T) are given by 

N 



D 



(b/e) 



1 

7ra 



n 



r 



(b/e 



■pN y(b/c)_ 



(24) 



2^1=1 



{b/e 



(25) 



where r[z] is the Gamma function. It thus follows that the 
TLL exponent associated with the i^-th band reads as 



N 



N 



1 



27rw.,- 



for the bulk position, and 

N 



1, (26) 



N 

j=i 



(27) 



for the edge. Equations ( |26] l and ( |27] i are the main findings of 
this article. 

Since D{u!,T, x) is given by the summation of the contribu- 
tions from each of the bands whose powers differ, the small- 
est values would be observed in actual experiments, such as 
photo-emissions and transport properties. 

Specifically, if TV = 2 with vpi ~ vp2 — vp, the present 
model is reduced to the conventional ID electron system, 
where the backward scattering between the different spins and 
the Umklapp scattering are ignored. In fact, by parameterizing 
as 



52(1,1) 


= .92(2,2) 


= 52|| -5i||, 


(28) 


ff2(l,2) 


= 52(2,1) 


= .92±, 


(29) 


54(1,1) 


= 54(2,2) 


= 54||, 


(30) 


54(1,2) 


= .94(2,1) 


= .94±, 


(31) 



Eqs. ( |26] l and jZTl ) lead to the familiar forms 



A('=)(l) = A(^)(2) = i(A7i+A7i)-l, 



(32) 
(33) 



with 



2ttvy - 


H54II 


+ 54± 


- .9211 - .92± + .91 II 


27ruF ~ 


H54II 


+ 54_L 


+ .9211 + 52± - .91 II 


2nvY - 


H54II 


- 54_L 


-52|| +52± +5i|| 


2'KV-p - 


\- .9411 


- .94± 


+ 52 II - .92± - 5i|| 



(34) 



(35) 



IV. W-CHANNEL SPINLESS FERMIONS 

In this section, we derive the matrix elements of the mutual 
interactions, which are included in the matrices K and V in 
Eqs. dU and (|9]l. As a simple example, we consider a quasi- 
ID spinless Fermion system where ID energy bands cross 
the Fermi energy, E-p . In addition, in order to clarify the ef- 
fects of the number of channels on the exponents, those for 
the case of long-range interaction are derived. 

The mutual interaction "Hint of the spinless Fermion can be 
expressed generally as 

Hint = ^ llxx'i'^{x)^^{x')V{\x-x'mx')iix), (36) 

with il.'{x) being the annihilation operator of the spinless 
Fermion. Since we are discussing low-energy physics, the 
interaction processes among the particles close to Ep are nec- 
essary. In order to obtain such interaction processes, the oper- 
ator ip{x) is expanded, using the eigenfunctions of the states 
across Ep, (t>v.K{x), as 



N 



Ipix) = E E '^'^,K4'iy,Kix), 
v=l K 



(37) 



where a^,K is the operator of the spinless Fermion with the 
eigenstate {v, K). By inserting Eq. (l37T i into Eq. ( |36] |, "Hint is 
expressed as 



\ E 



(38) 



where the matrix element of the mutual interaction is written 

as 



vi Kx,U2K2;i'3K3,ViKi 



- j I dxAyV{\x 



V\) 



X (t)t,.K,ix)K2,KAy)'t^'^3,KAy)(t>t^4.KAx)- 

(39) 

We note that as a result of momentum conservation, the re- 
lation Ki + K2 — K3 — K4 = nG holds, where G is the 
reciprocal lattice vector and n is an integer In the following, 
we discuss the case where the filling of each band is incom- 
mensurate. Then, only the normal processes satisfying n ~ Q 
are taken into account as 



UiKi,U2K2;U3K3,U4K4 — OK1+K2.K3+K4 



X K 



(K,,K2;K3,K4). 



(40) 
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In this case, 'Hint is expressed by 

^int =- 



1^1,1^2, 1^3, 1^4 Pi ,P2,P3,: 



'^1 i''2;f3,'^4 



^ "-Ii ,pi , i^l ,P2 , 1^2 '''=3 ,P3 , i^3 '^'=4 ,P4 , ) 



(41) 



where 

92{v, v') =LVuy-u'M{kFy, -kYu';-kYu',kYu), 
gi{v, v') =L{Vu,y-yM{kYv, k^v] k^u, kYu)5uy 

^v' \u' .uykFu ; kYif' 7 kpi/' , kpi/ ) 

- Viyy^uyikpiy, fepiy'; fcpiy, fcFi/')](l ' 



(46) 
(47) 



<5...')}- 
(48) 



where Cp,fe,^ = Oi^^p^p^ Here, the relations = p^fep^, + 
ki and jfeil ^ fcp^^ = I,-- - , A^) are used. Assuming 
kpi/ 7^ kpi^' for 7^ i^', Eq. (l4Tl i is written as "Hint = "Hint,! + 
Hint, 2 + Hint,4, where 



"Hi, 



1 



EE 



Here, the relation </>* f^{x) = 4)y-K{x), which is a result of 
time-reversal symmetry, is used. Note that gi (v, v'), 92 {v, v'), 
and g4^{v,v') are the real symmetric matrices. By compar- 
ing Eq. (l45T l with Eq. (IDi, we obtain 32(1', J^') = 52 (j^, i^') — 
and 54(1', !^') = .94(1', !^')- 
Here, we consider the case where the Fermi velocities of 
all the channels are equal to each other, i.e., upi/ = wp- In 
V^y-^y{pkp^,-pk-pv']pk-pv,-pkY,y'Y^^'^^'^on, we assume that the matrix elements are given by 

92{v, v') = gi{v, v') = g and ,91(1', v') = 0, which is a sim- 
ple but effective approximation for long-range interactions.— 
Here , the velocities of the excitation are obtained as vi = 

2,-. 



N 

k.k' ,q p—± u,^' — ! 

^ ^\+q,p,v'^\' -q.-py'^k' ,P,v'^f'^-py ^ 

N 



(42) 



Hint, 2 "2 E/ E/ E! ^t^,i^';i^',i^(pfcFiy, ~pkFv''i —pfcpt/' i pfcFiy)eigenvector corresponding to w = vi is written as Xi = 

k.k' ,qp=± v,v' = l 



?;f\/1 + Ng/ (irvp) and = (.7 
nd 



'"k-\-q,p,iy k' —q.—p^u' 
N 



(43) 



'Hint,4 =-| ^ y^Vv,y;u,y{pkFyiPkFiy]pkFu,pkFiy) 



k.k' ,q p=± v=l 

^ '^k+q,p,v^\' -q.p.v'^^' ,P,'^'^k,p,v 

N 



yj2nvF/N{l,--- ,1) 
tained as follows: 



N). The 

as J 

As a result, the exponents are ob- 



1 

1 

TV 



Wp 



Wp 



1 



(49) 
(50) 



i ^ ^ ^ |K,i.';i.',i.(pfcpi.,pfcp^/;pfcpi./,pfcpi.) Thus, the TLL exponents for both locations decrease with in- 
ky,qp=± v^v' creasing A^, and are proportional to A^~^/^ for iV ^ 1. 



' {pkpv , pkpu' ; pkpi' , P^Fi. 



^ ^\+q,p,v'^\' -q.py'^k' ,p,vCk^py^ . 



(44) 



Here, Hint.i represents the backward scattering, "Hint, 2 de- 
notes the forward scattering among the different branches, 
and Hint. 4 expresses the forward scattering between the same 
branches. It should be noted that we neglect accidental sit- 
uations in the momentum conservation, for example, fcp^^ — 
fcFi/2 = ~k-pu3 + fcp^4 with k-pui 7^ ^Fi/j and kp^^ ^ k-p^^, 
in the forward scattering among different branches. Equations 
(|42] |. (l43T i, and (gUi are reduced to 



1 



TV 



Hint E E E 

k.k' ,q p—:k — l 



[giiv. v')cl^gp ,^cl,_g _pyCk',pyCk-p,u 

+ 92{v, v')cl_^^^p ,^cl, _^ _pyCk' ,~pyCk,p,u 

+ 9i{v. i'')cl+g^p^A'-q,py'^k',pyCk,p,u}, (45) 



V. CONCLUSION 

In the present paper, we reformulated the TLL theory for 
multichannel ID Fermion systems. The theory obtained en- 
ables derivation of rigorous expressions for the local density 
of states and the corresponding TLL exponents, A^''^'^-' (:^), 
with respect to the i^-th band. The strategy for evaluating 
\{b/e) (^jy'j js summarized as follows: 

1. Define the functional forms of the ID eigenfunction 
<t>iy,K(x) and the interaction V{\x — y\) appropriate for 
the system being considered. 

2. Calculate KiXi, 1^2/^2; 1^3/^3, 1^4^:4 using Eq. 

3. Using the above result, set the mutual interaction terms, 
gi{v,v'), that are necessity to define the interaction 
Hint given by Eq. (|4]i. Particularly when consid- 
ering spinless Fermions, we can obtain gi{v,v') for 
i — 1,2,4 by substituting the results of step 2 into 
Eqs. 

4. Set (A'^^)y ,y/ and Vy,^i according to Eqs. (O and (|9]l. 
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5. Solve the eigenvalue problem ( fTsT l to obtain Vj and Xj 
foTj = !,■■■ ,N. 



6. Evaluate Y^ J and Y^^J from Eqs. (EB and (|22ll. 

7. Finally, we obtain the exponents X^''^i^) and X^'^^l^) 
from Eqs. (|26] | and ( l27l i. 

By applying the strategy for A^-channel spinless Fermions 
with long-range interactions, we have revealed that both TLL 
exponents approach zero in proportion to N~^^'^ for N ^ 
This finding implies that the power-law anomalies disap- 
pear and the energy-independent density of states, which is a 
manifestation of the realization of Fermi liquids, emerges for 
N oo. 

Before closing, we remark that the present theory began 
with the electronic Hamiltonian of the mutual interaction in 
Eq. which leads to the bosonic Hamiltonian Eq. (|7]i where 
no nonlinear terms exist. Even in the presence of interaction 
terms leading to nonlinear terms, the present theory can be 
useful if the the nonlinear terms are renormalized to zero. In 
such cases, it is necessary to take account of the renormaliza- 
tion of the parameters K and V by the diminishing nonlinear 
terms. 



with n being an arbitrary integer 

The mode expansion, together with the canonical quantiza- 
tion, leads to 

ej{x,t) = Cj + Q'j{x,t), (A5) 

(A6) 
(A7) 



e'j{x,t)^-l q 
^jix.t) = -- / q 



71" Jo ■' 

1 / COS qx 
■ 



where Cj is the c-number satisfying 
(l/x/2) Ef=i(^.,.7V^)Q = {n+ (1/2)}^, and b,{q) is 
the bosonic operator with [bj{q),b^j,{q')] = TTSjj'6{q — q'). 
The ultraviolet cutoff exp{~aq/2) is inserted if nec- 
essary in the g-integral in Eqs. iA6i and jAlj . Note 
that dtQj{x,t) = —Vjdx'^j{x,t). As a result of 
Eqs. ( IA5I )-( IA7I ). the Fermion field of the v-th band sat- 
isfies %lj^^^{x,t) ~ — x, t). The Hamiltonian is 
written as 
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Appendix A: Derivation of Eqs. fl9t-(l22l 

We discuss a semi-infinite system with its end at the origin. 
For convenience, we scale the bosonic fields as 



Qj{x,t) = ^jQj{x,t), 



(Al) 
(A2) 



where [Qj{x,t),^ji{y,t)] = \5jji9{x — y). By using field 
operators, the Hamiltonian is written as 

'H = jz'^ U{^ + {^.Q,f], (A3) 

where = —dx^j- The boundary condition at the origin 
requires the Fermion field for the i^-th subband i/j^ (0) = 0, 
i.e., ?/'-,i/(0) = —■)/'+, iy(0). This condition leads to 



3 = i 



(A8) 



The quantity 0), -0^ (x, i)}) is calculated as fol- 

lows: 



{{i'l{x,o),Mx,t)}) 

0), + -> -a;) 



— { l[G,^,{x,t)H,,^{x,t) + l[G,,,{x,t)H-j{x,t) \ 



27ra 

+ {x -x), 

where 



(A9) 



= exp <^(/^j(a;,0) - (x, t))^^ | , 

(AlO) 

H,,,{x,t) =exp|i [/,j(x,0),/,.,(a;,i)]| , (All) 



with = X^,,/^Q'.j{x,t) + 

2TT^j/vj{KX)i,,j^j{x,t). Here, we ignore the rapidly 
oscillating terms proportional to cxp(±i2fcFi/a;) because 
these contributions can probably not be observed directly 
due to averaging over several lattice sites in the experiments. 



7 



From Eqs. ( IA6I 1 and ( IA7l l. together with 

Gjy.j (a;, t) 



exp 



sm qx 



2 Jo 9 

X (1 - e"'"^**) (1 - e'"^"*) (1 + 2g{vjq)) 



2 7o 9 
X (1 - e--^«*) (1 - e'"^"^*) (1 + 2g{v,q)) 

Hy,j{x,t) 



(A12) 



exp <; ^ / q!l!L^ (ei-.9* _ c"'"^'?*) 



where 



(-'^t'j)^ 7-. UT^V\ l2 



= 2^w, (A14) 



and 5(e) — [(fl^ — 1) ^ is the Bose distribution function. As 
a result, 

({Vi(x,0),V.(2:,i)}> 
= — exp J - ^ [C+/,(0, t) + C_/,(a:, t)\ \ 



exp <! - ^ [C+ J, (0, t) + C_ J, (x, t)] 



(A15) 



where C± = {B^j ± A^j)/2 and 



(A13) 



■ 1 



2 log 



{vjt + 2xf 
{2xf 



log 



If sinh7rT(i + 2a;/?;j) 2t:Tx/z 
2 I 

+ I0, 



T(t + 2a;/i;j) sinh27rTx/uj 
sinh 7rr(t — 2xj Vj ) 2itTx / Vj 



Jj{x,t) 



TTT{t - 2x/vj) sinh27rra;/ 



■ 1 



Vi }' 



(A16) 



1 f, a - i(i;,i + 2a;) , a-i(vit-2x) 
= log ^^TT ^ + log ^- '- 

2 1 a — i2a; a + i2x 



(A17) 

Since Ij{oo,t) = Jj{oo,t) = 0, these results lead to 
Eqs. dull- 
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